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Abstract 

The phase transitions in a realistic system with triplet pairing, dense neutron 
matter, have been investigated. The spectrum of phases of the 3 i-2~ 3 -^2 model, 
which adequately describes pairing in this system, is analytically constructed 
with the aid of a separation method for solving BCS gap equations in states 
of arbitrary angular momentum. In addition to solutions involving a single 
value of the magnetic quantum number (and its negative), there exist ten real 
multicomponent solutions. Five of the corresponding angle-dependent order 
parameters have nodes, and five do not. In contrast to the case of superfluid 
3 He, transitions occur between phases with nodeless order parameters. The 
temperature dependence of the competition between the various phases is 
studied. 
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Superfluid systems manifesting triplet pairing between constituent spin- 1/2 fermions are 
exemplified by liquid 3 He at millikelvin temperatures and by the neutronic component of the 
quantum fluid interior of a neutron star. In the latter case, the density becomes high enough 
that the familiar singlet S'-wave gap has closed and pairing in the 3 P2 channel is favored by 
the spin-orbit force acting between neutrons The determination of superfluid phase 

diagrams for such systems over an extensive temperature range has presented a difficult 
challenge for theorists. On the one hand, the standard Ginzburg-Landau approach is valid 
only for temperatures near the critical temperature T c . On the other hand, iterative pro- 
cedures, commonly used to calculate the energy gap in systems with S'-wave pairing, suffer 
from slow convergence and uncertain accuracy when applied to the many coupled nonlinear 
integral equations that arise for BCS pairing in higher angular momentum states. 

Limitations of the standard iterative approaches are accentuated when one attempts to 
construct the superfluid phase diagram of the system, which is orchestrated by tiny energy 
splittings between the different solutions of the BCS pairing problem. This problem is 
ordinarily framed in terms of the set of equations Q 
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for the coupled partial-wave components that appear in the expansion 

Aa/3(p) = J2j,l,m ^l M (p) (^tj( n )) p °f the 2x2 gap matrix in terms of the spin-angle 
matrices (Gfj(n)) = Em s m l c l^ c iLM s M L Y LM L ( n )- The quasiparticle energy 
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is constructed from the gap components A J L M (p) together with the single-particle spec- 
trum £(p) ~ pf(p ~ Pf)/M* of the normal Fermi liquid, where M* is a suitable effec- 
tive mass. This quantity is rendered angle dependent by the spin trace S^ Jl 1 (n) = 



Tr 



G{ M (n)) Gi\ M Hu) 



further complicating explicit solution of the system (1). 



The pairing matrix elements {p\V^ L ,\pi) are generated by the spin-angle expansion 
V(p,Pi) = E i LVM(-l)^(p|^bi)G ! f J (n) (G^(nO)* of the block of Feynman dia- 
grams irreducible in the particle-particle channel. A salient feature of the vacuum nn in- 
teraction is that the components of the central forces nearly compensate each other, as 
evidenced in the experimental P-scattering phases ||. We shall assume that this feature 
is preserved by the effective nn interaction inside neutron matter. The dominant role of 
the spin-orbit component in promoting the 3 P 2 pairing channel at super-nuclear densities 
then implies that contributions to triplet pairing from nondiagonal terms with L', L\ 7^ 1 or 
J\ 7^ 2 on the r.h.s. of Eq. (1) can be evaluated within perturbation theory, in terms of the 
set of principal gap amplitudes A\ M (p), with M running from —2 to 2. In fact, time-reversal 
invariance implies that the problem may be treated in terms of only three complex functions, 
namely Af M (p) with M = 0,1,2. 

Essential further simplifications are possible if we apply the separation approach devel- 
oped in ||, which establishes the factorization 

Al M (p)=Dl M X (p), M = 0,l,2, (3) 

where xip) is a universal shape factor normalized to unity at p = pp. As a consequence, 
the elucidation of the phase diagram of dense superfluid neutron matter reduces to the 
determination of the three coefficients D\ M , since the character of the phase diagram itself 
is independent of the shape factor x(p)- 

Among the nondiagonal contributions to the r.h.s. of Eq. (1), two assume leading im- 
portance. The first contains the integral of the product V^S^ 2Ml Af 1 , while the second 
contains the integral of the product V^S^ 2Ml A% 1 . Considering only these contributions, 
we obtain the 3 P 2 - 3 F2 pairing problem, which has been treated numerically in earlier work, 
notably Refs. . The list of relevant states appears to be exhausted with the addition 

of the (less important) 3 Pq and 3 Pi pairing channels 0]. The rapid convergence of the non- 
diagonal integrals on the r.h.s. of the gap equations (1) greatly facilitates implementation 
of the perturbation strategy, since the overwhelming contributions to these integrals come 
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from the region adjacent to the Fermi surface. For E(jp) significantly in excess of the en- 
ergy gap value A F , the energies E(p) and practically coincide and, as a result, the 
angular integration yields zero. Thus, when dealing with the nondiagonal contributions, it 
is sufficient to know the minor gap components A\ M (p) (with M = 0,1,2) at the point 
p = pf, which may be efficiently evaluated in terms of the coefficients D\ M with the aid of 
the set ([!]) itself. In doing so, we retain on the r.h.s. of (P only the dominant contribution 
containing a large logarithmic factor ln(e^/ A^), where €f is the Fermi energy. This leads 
to the simple connection 



Al M ( P = p F )= V Dl M , (4) 

where 77 = — (pfI^iIIpf) A"f an d Vp = (pf\Vu\Pf) ■ Analogous linear relations hold for other 
minor components A^ M of the gap function (notably A^ and Aj M ). 

When calculated for in- vacuum neutron- neutron interactions 0, the r\ value depends 
smoothly on p, varying around 0.3 in the interval po < p < 3po [§]■ I n view of its relativis- 
ts origin, the spin-orbit force should not be much affected by polarization or correlation 
corrections, a judgment supported by the empirical analysis of spin-orbit splitting in finite 
nuclei. Medium modification of the tensor force may be more significant, especially in the 
environs of the pion- condensation phase transition [|K],[ll|]. Investigation of this possibility 



calls for special treatment and will be pursued elsewhere. 

A complete understanding of the phase diagram of systems that exhibit triplet pairing 
would require the inclusion of corrections to the in-medium interaction depending on the 
pairing gap itself. These "strong-coupling" corrections are most important near the critical 



temperature T c Hi2| , |13|| . Their effect is not addressed in the present work, which focuses on 
temperatures far enough below T c that pairing corrections to the vertex V (and hence to 
the parameter rj) may be ignored. 

We turn now to the determination of the key coefficients D\ M (M = 0, 1, 2) to leading 
perturbative order in the parameter r]. Making use of the connection (^), simple manipula- 
tions of Eqs. (1) at p = pf yield three coupled equations for these coefficients, 
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The search for solutions will be confined to those with real coefficients D\ since, as a 
rule, such states lie lower in energy than those with complex D\ M . Inserting the explicit 
form of S\i I2Ml (ja) into Eqs. (|5[) and @, we may arrive at a system of three equations 

A 2 + v F [A 2 (J + Js) ~ AiJi - J3] = ^f^2 , 
Ai + ujr [-(A 2 + 1) Ji + Ai(J + 4J 5 + 2J 3 )/4] = r/^n , 

l + w F [-(A 2 J 3 + A 1 Ji)/3 + J 5 ] =rrvFr , (7) 

for the two ratios Ai = D^/D^VE and A 2 = Df jDf\[E and the gap value A F . These 
equations are written so as to coincide at 77 = with the equations of the pure 3 P 2 pairing 
model solved in Ref. ||, where 

j* = y y /i(g,y)0(p) 2Eo(p) xfr)-^ (8) 

with /o = 1 - 3z 2 , f x = 3xz/2, f 3 = 3(2x 2 + z 2 - l)/2, and / 5 = (1 + 3z 2 )/2 and z = cosfl, 
a; = sin 9 cos 0, and y = sin 6 1 sin <fi. Of the integrals J, (with j = 1, • ■ ■ 5), only J5 contains a 
principal term going like ln(e F / A F ). 

Eqs. (0) have three familiar one-component solutions [|2],§|| corresponding to M = 0, 1, 
and 2. To establish the structure and the spectrum of the multicomponent solutions of the 
perturbed problem, we carry out a two-step transformation of the set (|7|). The integral J5 is 
responsible for introducing the gap value A F but is irrelevant to the phase structure. Thus, 
as a first step we combine the equations (0) so as to eliminate terms involving J5 from the 
first pair and, in addition, reduce the number of the Ji integrals in each equation to two. 
The resulting ^-independent equations are 

(A 2 + l)[3Ai(A 2 + 1) J - 2(A 2 - 2A 2 + 6) J x ] = r]Bx , 

(A 2 + 1)[(A 2 - 4A 2 ) Ji + Ai(A 2 + 1) J 3 ] = V B 2 , (9) 



where B 1 = 2Ai(2A 2 +3)r 2 -4(A 2 ; -3)ri-6Ai(A 2 +2)ro and B 2 = -Air 2 +4A 2 ri-3AiA 2 ro. The 
second step is taken under the assumption that A 2 ^ 1; otherwise the ensuing manipulations 
lose their meaning. As in Ref. || , we perform the rotation (x, z) = (— t sin d+u cos t cos i?+ 
Msin$), choosing the angle to remove the integral J\ from Eqs. (Q). In particular, if 
( = tan$ is chosen to obey the algebraic equation Ai£ 2 — (A 2 — 3)( — Ai = 0, substitution 
of the transformed integrals Jj into Eqs. @ yields 

(\ 2 + l)[A 1 J + A 2 J 3 ] = V B 1 , 

{\ 2 + l)[A 1 J + A 2 J 3 ] = -2r ] B 2 , (10) 

with A 1 = §A 1 (l + A 2 )(2-C 2 )-f(A 2 -2A 2 + 6)C and A 2 = -3A X (1 + A 2 )C 2 - (A 2 - 2A 2 , + 6)(. 

We observe that the left-hand members of these two equations (12) are identical. It is 
just this feature that leads to the universalities of the pure 3 P 2 pairing problem discovered 
in Ref. 0. Independently of temperature, density, and details of the in-medium interaction, 
the solutions of this restricted problem derived from Eqs. (9) at rj = (where the right-hand 
members of Eqs. (12) are trivially coincident) fall into two groups degenerate in energy, 
namely an upper group comprised of states whose angle- dependent order parameters have 
nodes and a lower group with nodeless order parameters (cf. Ref. |I4"|). In addition to the 



energy degeneracies, the multicomponent pairing solutions, which satisfy 

(A 2 + 2-2A 2 )(A 2 -2A 2 -6A 2 ) = 0, (11) 

display a degeneracy with respect to the coefficient ratios Ai and A 2 , since they generally 
define curves rather than points in the (Ai, A 2 ) plane. 

The strong parametric degeneracy inherent in pure 3 P 2 pairing is lifted in the case of 3 P 2 - 
3 F 2 pairing: at r] > the true solutions of the problem are represented by a set of isolated 
points in the (Ai, A 2 ) plane. Indeed, upon equating the right-hand members of Eqs. (|T0"|) in 
the small-// limit, one obtains an additional relation between the parameters \i(t] = 0) and 
A 2 (?7 = 0), viz. 

Air 2 - (A 2 - 3)n - 3Air = , (12) 



where the tm are defined by Eq. (§). This relation supplements the spectral condition (11) 
and lifts the parametric degeneracy. 

The system formed by ( |TT| ) and (|T2| ) is solved analytically by utilizing the same rotation 
in x — z coordinates as described above. After lengthy algebra, one arrives at the full 
set of solutions of the coupled-channel 3 P 2 - 3 F 2 pairing problem. A conspicuous feature of 
these solutions, made transparent by the separation method |J, is their virtually complete 
independence of the temperature T. 

Before cataloging and classifying the 3 P 2 - 3 F 2 pairing solutions, it is worth pointing out 
that the particular solution A2 = — 1 found in the pure 3 P 2 problem remains intact upon 
switching on the 3 F 2 channel, but it completely disappears if any further channel (e.g. 3 Po) 
enters the picture. 

Fig. 1 summarizes the results of our analytical considerations. Only the right half of 
the (Ai,A2) plane is shown, since the relevant energies are independent of the sign of Aj. 
Besides the well-known single-component solutions with \M\ = 0, 1, or 2, the collection 
of unitary solutions of the 3 P 2 ~ 3 F 2 pairing problem contains ten multicomponent solutions, 
corresponding to more complicated superfluid phases. Five of these additional solutions, 
denoted Ok (k = 1, . . . , 5), have nodeless order parameters and include: 

(a) Two two-component solutions 0±3, identical to those found in the pure 3 P 2 pairing 
problem, with Ai = and A 2 = ±3. 

(b) Three three-component solutions, two of which, Oi and O4, are associated with the 
upper branch of Xf — 2A 2 — 6A 2 = and have A 2 = 3(^/21— 4)/5 and A 2 = 3, respectively; 
while a third, 2 , is associated with the lower branch of the same equation and has 
A 2 = -3(V2l + 4)/5. 

The remaining five solutions of the problem, denoted by {k — 1, . . . , 5), are distinguished 
by order parameters that do have zeros. This set includes: 

(a) Two two-component solutions X±i, again identical to those found in the pure 3 P 2 
pairing problem, with Ai = and A 2 = ±1. 
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(b) Three three-component solutions, X 2 , X3, and X±, associated with the parabola A2 = 
Xl/2 + 1 and having A 2 = 13 — 2y / 35, A 2 = 3, and A 2 = 13 + 2\/35, respectively. 

These general features of the spectrum of solutions of the 3 P 2 - 3 F 2 problem are expected to 
persist even if rj is not so small. 

To construct the phase diagram of superfluid neutron matter, the third of Eqs. (|7|) must 
be brought into play to determine the gap values Ap for the different solutions. We must 
then compare the free-energy shifts F s = — J A 2 F (g)g~ 2 dg due to pairing in the corresponding 
superfluid states, where g is the relevant pairing coupling constant. 

At low T, the only true contestants for existence and ascendancy are the solutions with 
nodeless order parameters, since the other solutions lie too high in energy; the gap between 
the states of the two different groups cannot be bridged if the value of \r)\ stays rather 
small. Thus, even though the degeneracy of the 3 P 2 pairing problem in the Ai — A 2 plane, 
embodied in the relation (11), is entirely removed when the 3 P 2 - 3 F 2 coupling is switched on, 
the degeneracies in the energetic spectrum of the different superfluid phases are only partially 
lifted. Instead, this spectrum decays into several groups of nearly degenerate states, the 0\ 
and 2 phases forming the lowest-energy group, the next higher group being composed of 
the phases 0±% along with the one-component phase with M = 0, and so on. As seen in 
Fig. 2, the splitting between the two groups lowest in energy shrinks as T increases, until, at 
T ~ 0.7T C , their roles are interchanged and transitions occur. However, unlike the A-B phase 
transition in superfluid 3 He, where the order parameter of the A-phase solution has nodes 
while that of the B-phase does not, the transition in neutron matter takes place between 
phases with nodeless order parameters. Introduction of the 3 Po and 3 Pi pairing channels 
produces a further lifting of energy degeneracies and hence a further complication of the 
phase diagram of superfluid neutron matter. Unfortunately, it is not possible to resolve the 
resulting fine structure of the phase spectrum without more detailed information on the 
in-medium effective interaction between neutrons, which is currently unavailable. 

The temperature region near T c provides another venue for possible superfluid phase 



S 



transitions. In this regime, strong-coupling corrections should no longer be ignored ||13|| , and 
the character of the phase diagram may be influenced significantly by external magnetic 



fields. Finally, fermion condensation, occurring as a precursor to pion condensation [15 



should be mentioned as a source of exotic phase transitions in superfluid neuron matter at 
densities approaching the critical value for collapse of the pion mode JTU| . 

If, as we have shown, the phase diagram of dense neutron matter exhibits several triplet 
superfluid phases, then phase transitions between different phases are expected to occur as 
the neutron star cools. Since the gap value changes in these transitions, their occurrence 
may ultimately be detected in the thermal history and/or the rotational dynamics of the 
star, for example in variation of its moment of inertia or in alterations of the distribution of 
the angular momentum between the crust and the vortex system, resulting in a change of the 
star's angular velocity. Possible scenarios leading to distinctive observable manifestations 
will be presented separately. 
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Figure Captions 



Fig. 1. Multicomponent solutions of the 3 P 2 - 3 F 2 pairing problem. Solutions whose order 
parameter is nodeless [exhibit nodes] in the case of pure 3 P 2 pairing are indicated by filled 
[open] circles. 

Fig. 2. Difference 5A 2 F (T) between the A F (T) values for the phases 0\ or 2 and for the 
phases 0±s (or the one-component phase with M = 0), measured in units of \5A F (T = 0)|. 
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